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Abstract. The mass-loss rates and dust-to-gas ratios of 
obscured Asymptotic Giant Branch (AGB) stars are in- 
vestigated for samples with different initial metallicities: 
in the Small and Large Magellanic Clouds (SMC & LMC) 
and in the Milky Way. The properties of their circumstel- 
lar envelopes can be explained in a consistent way if, both 
for obscured M-type AGB stars and for obscured carbon 
stars, the total (gas+dust) mass-loss rate M depends only 
weakly on initial metallicity whilst the dust-to-gas ratio tp 
depends approximately linearly on initial metallicity. 

Key words: circumstellar matter - Stars: mass loss - 
Stars: AGB and post- AGB - Stars: winds, outflows - Mag- 
ellanic Clouds - Infrared: stars 



1. Introduction 

Asymptotic Giant Branch (AGB) stars develop strong 
mass loss at rates up to M ~ 10~ 5 M Q yr _1 or more. The 
optical light of AGB stars with the highest M is almost 
entirely absorbed by their dusty circumstellar envelopes 
(CSEs), and re-emitted at longer wavelengths. These ob- 
scured AGB stars become very bright infrared (IR) ob- 
jects, outshining any other star in a galaxy except for a 
few red supergiants (RSGs) . During a brief period in their 
lives, they loose 50 to 80% of their initial mass. This makes 
them important contributors to the chemical enrichment 
of the interstellar medium (ISM) — they are possibly the 
main sources of dust particles in the Universe. 

Individual galaxies differ in their metal content as a 
result of different star formation histories. For instance, 
two of our nearest neighbours, the Large and Small Mag- 
ellanic Clouds (LMC & SMC) have current metallicities 
a factor ~ 2 and ~ 5 less than the average metallicity in 
the Milky Way (~solar). When pursuing a quantitative 
description of the history of star formation and chemical 
enrichment within a galaxy, it is essential to correctly take 



into account the mass loss from AGB stars and RSGs. This 
requires the dust-to-gas ratios in their CSEs to be known. 

The paper is organised as follows: in Sect. 2 formulae 
are derived for deriving mass-loss rates and dust-to-gas ra- 
tios from measurements of optical depths and either lumi- 
nosities or expansion velocities. Sect. 3 introduces samples 
of obscured AGB stars in the MCs and in the Milky Way 
and addresses the near-IR colours of their circumstellar 
envelopes. Relative mass-loss rates and dust-to-gas ratios 
are determined in Sect. 4 for the magellanic and galac- 
tic circumstellar envelopes around obscured AGB stars, 
and the results are discussed in Sect. 5. In Appendix A 
new identifications of mass-losing AGB stars with IRAS 
point sources in the LMC are presented and in Appendix 
B expansion velocities are discussed. 

2. Formulae for dust-driven winds 

The optical depth r is proportional to the column density 
and dust-to-gas mass ratio ip- The column density is pro- 
portional to the total mass density p and inner radius Ri 
of the dusty CSE. The dust grains in the dust-formation 
zone are assumed to be in radiative equilibrium with the 
incident stellar radiation field, keeping the effective tem- 
perature T e ff and the dust condensation temperature 
fixed. Then Rf oc L, the stellar luminosity. The continuity 
equation yields M = 47rr 2 /9i> exp , with gas+dust mass-loss 
rate M, and expansion velocity v exp . Here the drift ve- 
locity — the velocity difference between the gas and dust 
fluids — is neglected, i.e. gas and dust are assumed to be 
well coupled (cf. Lamers & Cassinclli 1999). Gas pressure 
and wind-driving mechanisms other than radiation pres- 
sure on dust are neglected as well, although they become 
important for wind speeds below a few km s _1 (e.g. Stcffcn 
et al. 1997, 1998). It follows that the optical depth 



r(A) cx 



u cx.p 



(i) 



Send offprint requests to: Jacco van Loon, jacco@ast.cam.ac.uk 
* This paper is based in part on data obtained at the South- 
African Astronomical Observatory. 



The constant of proportionality includes a factor k(A), the 
wavelength-dependent opacity of the dust, and a temper- 
ature dependence as (T d /T e g) 2 . 
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In a radiation-driven outflow the matter-momentum 
flux is coupled with the stellar photon-momentum flux 
via the momentum equation (Gail & Sedlmayr 1986) 



Mi> exp cx tL 



(2) 



with t the flux-weighted optical depth. In general, the 
ratio of t and r(A) depends on the mass-loss rate. De- 
tailed computations such as those presented in Habing et 
al. (1994), however, show that this dependence vanishes 
for M > M yr -1 , which applies to the class of 

obscured AGB stars under study here (van Loon et al. 
1999b). Combination of Eqs. (1) and (2) then leads to a 
description of the expansion velocity in terms of dust-to- 
gas ratio and luminosity: 



(3) 



Eq. (2) can also be used to eliminate w cxp from Eq. (1), 
yielding a relation between the mass-loss rate and dust- 
to-gas ratio, and the optical depth and luminosity: 



logM + 0.5 log tp + constant = logr + 0.75 log L 



(4) 



Alternatively, Eq. (2) can be used to eliminate L from Eq. 
(1), yielding a relation between the mass-loss rate and 
dust-to-gas ratio, and the optical depth and expansion 
velocity: 



log M + 2 log V> + constant = log r + 3 log v c 



(5) 



For magellanic stars luminosities are easier to measure 
than expansion velocities, and for these stars it is advan- 
tageous to make use of Eq. (4) . For galactic stars the op- 
posite is true, and for them Eq. (5) is the formula to use. 

The constants in Eqs. (4) and (5) are related. Calling 
the constants of proportionality of Eqs. (1) and (2) respec- 
tively a and 0, the constant in Eq. (4) equals \og{a/0) and 
the constant in Eq. (5) equals \og(a 2 0). The constant of 
proportionality in Eq. (3) equals \fa0. The values of a 
and depend on the properties of the dust species, and 
calibrating them is an important yet extremely difficult 
task that will not be exercised here. 



3. Samples of obscured AGB stars 

The optical depth scales with near-IR colours: r cx (H — 
K) — (H — K)q. The stellar photosphere is assumed to 
have (H — K) = 0.2 mag. Fig. 3 in paper IV shows that 
galactic M-stars have (H — K) ~ 0.2 mag, whilst galactic 
carbon stars may be slightly redder with (H — K)q ~ 0.4 
mag. The sensitivity to the choice of (H — K) rapidly 
vanishes as stars become obscured at (H — K)^l mag. 

Fig. 1 shows the (H — K) colours of the obscured AGB 
stars and RSGs in the MCs (Loup et al. 1997; Zijlstra et al. 
1996; van Loon et al. 1997, 1998a: papers I to IV). These 
samples are based on the identification of optical and near- 
IR counterparts of point sources detected at 12 and/or 25 
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Fig. 1. (H — K) colour versus bolometric magnitude for 
obscured AGB stars and red supcrgiants in the Magellanic 
Clouds. Colours are bluer for the more luminous stars. 



/im by IRAS. The SMC photometry is mainly from Groe- 
newegen & Blommaert (1998). Bolometric magnitudes for 
the SMC stars were determined in the same way as for 
the LMC stars by spline fitting to the spectro-photometric 
energy distribution (see Whitelock et al. 1994), adopting 
distance moduli of 18.55 and 18.97 mag for the LMC and 
SMC, respectively (Walker 1999). Also included are two 
newly identified obscured AGB stars in the LMC (Ap- 
pendix A). 

The reddest, optically thickest stars, with (H — K)^3 
mag, are found among AGB stars with Mboi ~ —5 to —6 
mag. No such red objects are known in the MCs among the 
brightest AGB stars with Mboi ~ — 7 mag, nor amongst 
the RSGs. This is partly because bolometrically fainter 
stars have smaller inner radii of the dusty CSE and smaller 
expansion velocities, yielding larger t (Eq. (1)) and redder 
(H — K) colours at a given mass-loss rate. 

Whitelock et al. (1994) have searched for Long Period 
Variables (LPVs) in the South Galactic Cap (SGC). Their 
sample consists of both optically bright and obscured AGB 
stars. Their (H — K) colours are plotted versus bolometric 
luminosity in Fig. 2 (the Mira P-L relation was applied), 
where different symbols are used according to the 25 /im 
flux density measured if the star were at the distance of 
the LMC. The obscured stars that are detected by IRAS 
in the LMC have 525^0.1 Jy. Their SGC equivalents have 
very red (H — K) colours, though not redder than the 
reddest in the LMC (crosses). Optically bright AGB stars 
in the SGC sample have typically (H — K) ~ 0.5 and 
Mboi ~ —4 to —5 mag. The MC samples do not contain 
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Fig. 2. (H — K) colour versus bolometric magnitude for 
AGB stars in the South Galactic Cap (Whitclock ct al. 
1994). Different symbols are used according to the 25 /zm 
flux density measured if the star were at the distance of 
the LMC. Crosses represent obscured AGB stars in the 
LMC (from Fig. 1). 



such objects because their mid-IR emission is too faint 
to have been detected by IRAS at the distances of the 
MCs. On the other hand, the SGC sample is devoid of 
the brightest AGB stars with Mboi ~ — 7 mag as well as 
RSGs, because in the Milky Way such massive stars are 
preferentially found in the galactic plane. 

There is a cluster of SGC stars with (H - K) ~ 0.8 
and Mboi ~ —5.2 mag (Fig. 2). These stars clearly show 
circumstellar reddening, but the mid-IR emission from 
their CSEs is just below the detection limit of IRAS when 
placed at the distance of the LMC. This leaves open the 
possibility of the existence in the LMC of a potentially 
large population of AGB stars with moderate mass-loss 
rates and luminosities. Indeed, in paper III several field 
stars were found that are not related to a nearby IRAS 
source but that nevertheless had near-IR colours indica- 
tive of reddening. Recent ISO observations confirm the 
presence of this AGB population (Loup et al. 1999). 

Wood et al. (1998) find LPVs in the Galactic Centre 
with K-band magnitudes from 5 to 13 after correction for 
interstellar extinction. At the distance of the LMC this 
would yield K-band magnitudes from 9 to 17, i.e. within 
the sensitivity of the searches in papers II & III and in 
Groenewegen & Blommaert (1998). The interstellar ex- 
tinction corrected (H — K) colours of the Galactic Centre 
LPVs average ~ 2 mag and are ~ 4 mag maximum, similar 
to the (H — K) colours of the obscured AGB stars in the 



LMC. Wood et al. argue that their sample includes stars 
with initial mctallicities a few times solar, and Blommaert 
et al. (1998) indeed find very red objects with (H — K) > 4 
mag as inferred from their K- and L-band photometry. 

Groenewegen et al. (1998) observed and modelled ob- 
scured carbon stars in the Milky Way. They used the P-L 
relation for carbon Miras to infer distances to the indi- 
vidual stars, that are found to be typically within 2 kpc 
from the Sun. Their two most obscured carbon stars have 
(H — K) = 5.5 and 8.0 mag. They also compiled a sample 
of oxygen-rich stars with near-IR photometry, pulsation 
periods and expansion velocities, without overlap with the 
SGC and Galactic Centre samples. The most obscured of 
these M stars have (H — K) ~ 6 mag. 

4. Relative mass-loss rates and dust-to-gas ratios 

In determining mass-loss rates and dust-to-gas ratios in 
the MCs one can rely on the use of luminosities and 
eliminate v cxp from Eq. (1). Luminosities are difficult to 
measure for stars in the Milky Way, however, due to un- 
known distances and severe interstellar extinction. Often 
the Mira P-L relation is applied to infer distances to in- 
dividual stars, but this relation possibly breaks down for 
stars whose stellar mantles have been significantly reduced 
due to mass loss (e.g. Blommaert ct al. 1998; Wood et al. 
1998; Wood 1998), or obscured AGB stars may fall on 
other, parallel sequences (Bedding & Zijlstra 1998; Wood 
1999). For the Milky Way stars, L is eliminated from Eq. 
(1), leaving u exp to be measured. Because v exp has been 
measured for only a few LMC stars, Eq. (3) will be used 
for the remaining LMC stars to estimate v exp from L (see 
Appendix B). Although these estimates are likely to be 
accurate within a few km s _1 , the contribution of the un- 
certainty in v cxp is enlarged by a factor three in Eq. (5). 

Eqs. (4) & (5) are used here to compare the combina- 
tion of mass-loss rate and dust-to-gas ratio between the 
MCs and between the LMC and the Milky Way, respec- 
tively. If the mass-loss rate depends on luminosity but not 
on initial metallicity and if the luminosity distributions of 
the stars in these samples are identical — implying identi- 
cal star formation histories — then the distributions over 
the combination of M and ip are identical except for pos- 
sible offsets due to different mean values for ip among the 
different samples. The stars within each of the samples 
are likely to cover a range in initial mctallicities due to 
their different progenitor masses and hence different for- 
mation epochs. Similar star formation histories, however, 
are anticipated to result in similar distributions over ini- 
tial metallicity, and it is therefore meaningful to assign a 
mean initial metallicity and a mean value for tp to each of 
the samples. 

Eq. (4) is used to calculate the combination of mass- 
loss rate and dust-to-gas ratio from the optical depths and 
luminosities for the magellanic stars (Fig. 3). The stars 
with Mb i < —7.5 and > —7.5 mag are considered to be 
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Table 1. Mean and error in the mean of the values of (log M + 0.5 log ip ± constant) for the distributions of obscured 
AGB stars in the LMC and SMC. Also the differences between the various (sub-)distributions are listed. 





[total] 


[C] 


[M] 


[C] - [total] 


[M] - [total] 


[C] - [M] 


[LMC] 


-0.85 ±0.04 


-0.83 ±0.05 


-0.83 ±0.09 


0.02 ± 0.06 


0.02 ± 0.09 


0.00 ±0.10 


[SMC] 


-1.04 ±0.05 


-1.14 ±0.16 


-1.10 ±0.36 


-0.10 ±0.17 


-0.06 ±0.36 


-0.04 ±0.39 


[LMC] - [SMC] 


0.19 ±0.07 


0.32 ±0.17 


0.28 ±0.37 
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Fig. 3. Mass-loss rate M and dust-to-gas ratio i/j calcu- 
lated from Eq. (4), versus Mboi- There is a maximum to 
M, increasing with higher luminosity. Mass-loss rates may 
be the same in SMC and LMC if tp is sufficiently lower in 
the SMC. 



RSGs and obscured AGB stars, respectively (see van Loon 
et al. 1999b), except for the well-studied luminous ob- 
scured AGB star IRAS05298-6957 (M bo] = -7.65 mag) 
and the weakly mass-losing RSG HV2700 (M bo i = -7.45 
mag). IRAS04498-6842 (M bo i = -8.28 mag) was classi- 
fied as an AGB star (paper II), but here it is re-classified 
as a RSG (see also paper IV). Fig. 4 shows the cumulative 
distributions (normalised to unity) of the obscured AGB 
stars in the LMC (solid) and SMC (dotted) over the value 
of (logM±0.51ogV'±constant). These do not include the 
few stars that have lower limits to their (H — K) colours. 
The distributions for the subsamples of spectroscopically 
confirmed carbon stars and oxygen-rich (M-type) stars 
(van Loon et al. 1999b, and references therein; Groenewe- 
gen & Blommaert 1998) are boldfaced in Figs. 4a and 
b, respectively. Despite the small numbers for some of the 
sub-distributions, especially the M-type stars in the SMC, 
the shapes of the cumulative distributions are very sim- 
ilar. Both in the LMC and SMC the obscured M-type 
AGB stars and the obscured carbon stars have distribu- 




Fig. 4. Normalised distribution of obscured AGB stars 
(MboiKi — 7.5 mag) over mass-loss rate M and dust-to-gas 
ratio ip. Boldfaced are the sub-distributions of confirmed 
M-type stars (a) and confirmed carbon stars (b) amongst 
the obscured AGB stars. 



tions that are coincident to a very high degree (Table 1). 
The only difference in the distributions is the systematic 
offset between those in the LMC compared to those in the 
SMC (Table 1). This offset is a 3(7 result for the obscured 
AGB stars as a whole, a 2a result for the obscured carbon 
stars, and consistent for the obscured M-type AGB stars 
(though not significant by itself). 

Eq. (5) is used to calculate the combination of mass- 
loss rate and dust-to-gas ratio from the optical depths 
and expansion velocities for the obscured AGB stars in 
the LMC and the Milky Way (Fig. 5). Only LMC stars 
for which periods are known from the literature (Wood 
et al. 1992; Wood 1998) are plotted. Stars in the differ- 
ent samples are not similarly distributed over P, with the 
LMC sample having an average pulsation period longer 
than that of the Milky Way sample, possibly due to selec- 
tion effects or differences in the star formation histories. 
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Fig. 5. Mass-loss rate M and dust-to-gas ratio tp calcu- 
lated from Eq. (5), plotted versus pulsation period P for 
obscured AGB stars in the LMC (Wood et al. 1992; Wood 
1998; van Loon ct al. 1998b), South Galactic Cap (White- 
lock et al. 1994), Galactic Centre (Wood et al. 1998) and 
"Solar neighbourhood" (Groenewegen et al. 1998). Dot- 
ted vertical lines delineate regimes of periods P < 500, 
500<P < 800 and P > 800 d. 

Hence the populations of stars in these samples may not 
be directly compared. 

Two regimes of pulsation period will be focussed on: 
one for 500<P < 800 d, and another for P > 800 d. 
The first group consists of stars that evolve beyond the 
optically bright Mira phase with moderate M and into 
the obscured AGB phase (Jura 1986). The distribution of 
these stars in Fig. 5 suggests that they still increase in 
mass- loss rate and/or that they represent a large range in 
stellar masses. The second group consists of more evolved 
(or more massive) stars for which M has become so high 
that they shed their mantles on a timescale much shorter 
than the nuclear burning timescale. Hence these stars stay 
at constant luminosity while their pulsation periods keep 
increasing as their mantles are steadily diminished. The 
lack of any clear correlation between M and P in this 
part of Fig. 5 suggests that thereby the mass-loss rate 
remains essentially constant. Stars in the LMC sample are 
predominantly found in the second class of objects, whilst 
these kind of stars are very rare in the Galactic Centre 
sample due to their faintness in the K-band. M-type stars 
in the "Solar neighbourhood" sample of Groenewegen et 
al. (1998) populate both regimes rather evenly. 

The relative distributions of mass-loss rates and dust- 
to-gas ratios are derived from Fig. 6 for the LMC (dot- 
ted), "Solar neighbourhood" (M-type; solid) and Galactic 
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Fig. 6. Normalised distribution of obscured AGB stars in 
the LMC and in the Milky Way over mass-loss rate M 
and dust-to-gas ratio tp for pulsation periods P > 800 d 
(a) and 500<P < 800 d (b). 

Table 2. Mean and error in the mean of the values of 
(log M ± 2 log ip ± constant) for the distributions of ob- 
scured AGB stars in the LMC, in the Galactic Centre 
(GC), and in the "Solar neighbourhood" (Groen; exclud- 
ing the point at P = 802 d and a value of —2.7). Also 
the differences between the various (sub-)distributions are 
listed. 



[C] 



[M] 



[C] - [M] 



500<P < 800 d 
-1.16 ±0.04 

-0.06 ±0.10 -0.55 ± 0.08 0.49 ±0.13 
-0.19 ±0.08 

-[LMC] 1.10 ±0.11 

[Groen] 0.37 ±0.12 



[LMC] 

[Groen 

[GC] 

[Groen 

[GC]- 



[LMC] 

[Groen] 

[GC] 

[Groen] 

[GC] - 

[GC]- 



- [LMC] 

[LMC] 

[Groen] 



P > 800 d 
-0.75 ±0.05 -0.79 ±0.12 
0.17 ± 0.24 0.01 ±0.06 
0.15 ±0.06 
0.93 ± 0.25 0.79 ±0.14 
0.93 ±0.14 
0.14 ±0.09 



0.04 ±0.13 
0.17 ±0.25 



Centre (dashed), in the same way as before by estimat- 
ing the relative offsets in (log M + 2 log ip ± constant) of 
the normalised cumulative distributions. The sample of 
obscured AGB stars with 500<P < 800 d in the "So- 
lar neighbourhood" (Groenewegen et al. 1998) is clearly 
inhomogencous: (i) the distribution of obscured M-type 
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AGB stars is much natter than all other distributions, in- 
dicating a wide range of stellar parameters, and (ii) the 
distribution of obscured carbon stars is offset with respect 
to the distribution of obscured M-type AGB stars — a 4<r 
result (Table 2). This may reflect differences in age (initial 
stellar mass) and/or initial metallicity, as the sub-samples 
of obscured M-type AGB stars and obscured carbon stars 
in the Groenewegen et al. sample have been selected in 
a very different way. The distributions of the obscured 
M-type AGB stars and obscured carbon stars in the LMC 
sample are, again, indistinguishable (Table 2). The Galac- 
tic Centre sample is positively offset with respect to both 
the Groenewegen et al. sample (3cr result for 500<P < 800 
d) and the LMC sample (7a result). 



5. Discussion 

5.1. Metallicities of obscured AGB stars in the 
Magellanic Clouds and in the Milky Way 

What are the initial metallicities of the obscured AGB 
stars in the samples under study? First I discuss the dust- 
to-gas ratios in the ISM, before discussing metallicities of 
relatively young stellar populations. These considerations 
lead to estimates of the typical initial metallicities of the 
obscured AGB stars under study. These will then be used 
to correlate the mass-loss rates and dust-to-gas ratios with 
initial metallicity. 

Most of the mass of stars that eventually evolve into 
AGB stars is returned to the ISM, enriched with dust. 
These intermediate-mass stars are sufficiently numerous 
and short-lived that the ISM has been recycled by them at 
least several times over the history of their parent galaxy. 
A 5 M Q star has a lifetime of <~ 1.1 x 10 8 yr (Marigo et 
al. 1999), and there have been already <~ 10 2 generations 
of these massive AGB stars at work in mixing dust with 
the ISM. Dust-to-gas ratios in the ISM of the LMC are 
found to be ~ | of that of the ISM in the Solar neigh- 
bourhood (van Gendercn 1970; Koornneef 1982; Clayton 
& Martin 1985), and in the SMC it is ~ ^ of that in the 
Solar neighbourhood (van den Bergh 1968; van Genderen 
1970; Lequeux et al. 1982; Bouchet et al. 1985). Dust is 
being destroyed in the ISM by shocks and gas accretion, 
especially in Star Formation Regions. Hence the dust-to- 
gas ratios in the ISM may only pose a lower limit to the 
dust-to-gas ratios in the CSEs of obscured AGB stars. 

Initial metallicities for relatively young (10 7 ~ 8 yr) field 
stars are found to be about 2 and 5 times lower in the 
LMC and SMC, respectively, compared to Solar metallic- 
ity (Spite et al. 1989a,b; Russell & Bessell 1989; Meliani et 
al. 1995; Luck et al. 1998). When depicting the chemical 
evolution of the MCs, however, it is clear that the initial 
metallicity with which stars were born a few Gyr ago was 
considerably lower — roughly twice — than what is mea- 
sured in massive stars today (de Freitas Pacheco et al. 
1998; Da Costa & Hatzidimitriou 1998; Bica et al. 1998). 
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Fig. 7. Mass-loss rates and dust-to-gas ratios for the ob- 
scured AGB stars in the SMC, LMC, "Solar neighbour- 
hood" and Galactic Centre, as a function of their initial 
metallicities. 



Obscured AGB stars, with a characteristic age of 10 9 yr, 
are thus expected to have initial metallicities somewhat 
lower than those measured in young field stars. Stars in 
the central regions of the Milky Way galaxy are found to 
cover a range in metallicity from sub- to super-solar, with 
the most easily observed obscured AGB population likely 
to comprise relatively massive stars of slightly super-solar 
initial metallicity (Rich 1988; Wood et al. 1998). 

These considerations lead us to adopt a typical initial 
metallicity (logarithmic) of [Z/H] = -0.85 ± 0.15 for ob- 
scured AGB stars in the SMC, [Z/H] = -0.45 ± 0.15 in 
the LMC, [Z/H] = 0.00±0.10 in the Solar neighbourhood, 
and [Z/H] = 0.20 ±0.20 in the Galactic Centre. The mar- 
gins are rough estimates for the typical range in initial 
metallicities, meaning that the typical metallicity for the 
sample is not expected to lie outside of these margins. 

5.2. Mass-loss rates and dust-to-gas ratios 

The values for the combination of mass-loss rates and 
dust-to-gas ratios as derived from using Eqs. (4) & (5) 
and listed in Tables 1 & 2 are plotted in Fig. 7. First or- 
der polynomials are drawn through the points. For the 
obscured M-type stars with 500 <P < 800 d the polyno- 
mial was forced to cross the LMC data point, while the 
average was taken of the slopes obtained by either omit- 
ting the "Solar neighbourhood" data point or the Galactic 
Centre data point. The slopes of the polynomials can be 
used to constrain the dependence of the mass-loss rate 
on initial metallicity z (with logz = [Z/H]) and the de- 
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pcndencc of the dust-to-gas ratio on z, assuming that 
M oc z a and oc z b . It then follows that the slope for 
(logM + 0.51ogV' + constant) equals a+0.56, and the slope 
for (logM + 21og-0 + constant) equals a + 2b. The sub- 
samples that had been selected according to the pulsation 
periods have, in fact, very similar slopes, and hence their 
averages are taken. Thus I find for the obscured M-type 
AGB stars that a + 0.5b ~ 0.68 and a + 2b ~ 1.71, whilst 
I find for the obscured carbon stars that a + 0.56 ~ 0.78 
and a + 2b ~ 2.24. This yields for the dependences of the 
mass-loss rate and dust-to-gas ratio on initial metallicity 
for obscured AGB stars: 

log Mm = 0-3li.3 1°§ Z + constant 

logMc = log z + constant 

logV'M = O-Tto'J logz + constant 

log^c = 1-0^0.3 1°§ Z + constant 

The reader should realise that these exponents are indica- 
tive, but not very accurate. The errors (excluding those 
resulting from uncertainties in the adopted initial metal- 
licitics) mainly result from the small number of spectro- 
scopically confirmed carbon and M-type stars in the SMC. 
If the initial metallicity dependence of the mass-loss rate 
and dust-to-gas ratio is really as similar as suggested here, 
then the obscured AGB stars in the SMC and LMC may 
be compared without distinction by chemical type, and 
the resulting errors on the exponents will become much 
smaller without significantly changing the values of the ex- 
ponents themselves. Clearly, it is shown that the method 
employed here has the prospect of constraining the initial 
metallicity dependences of the mass-loss rate and dust-to- 
gas ratio for obscured AGB stars. Suitable data have only 
recently become available, and more such data are needed 
to improve on the preliminary results. 

5.3. Obscured M-type AGB stars and carbon stars 

Obscured M-type AGB stars and obscured carbon stars 
have been treated separately for two reasons: (i) the op- 
tical properties of their circumstellar dust are different, 
possibly leading to differences in the constants of propor- 
tionality a and [3 in Eqs. (1) and (2), respectively, and (ii) 
their dust-to-gas ratios and/or mass-loss rates may de- 
pend differently on initial metallicity. I here discuss these 
two issues in more detail. 

Both a and j3 depend on the dust-type through the 
wavelength dependent opacity k(X) and the flux-weighted 
opacity k, respectively. Using Eq. (4) in the LMC yields 
indistinguishable distributions of the obscured M-type 
AGB stars and of the obscured carbon stars. However, 
in the LMC, obscured carbon stars are bolometrically 
fainter and exhibit lower mass-loss rates than obscured 
M-type AGB stars (van Loon et al. 1999b), which is ex- 
pected to result in an offset between the distributions over 
(logM + 0.5 log V> + constant). The fact that this offset is 



not seen must mean that, by mere coincidence, the dif- 
ferences in (logM -I- 0.5 log VO are counteracted upon by 
the differences in k(X)/k: at a given luminosity, dust-to- 
gas ratio and mass-loss rate, obscured carbon stars are 
redder than obscured M-type AGB stars. Also the coinci- 
dence between the distributions of obscured M-type AGB 
stars and obscured carbon stars using Eq. (5) in the LMC 
for 500<P < 800 d must be coincidental and/or due to 
low number statistics (only two obscured carbon stars), 
as the constant in Eq. (5) is proportional to k 2 (X)k and 
differences in optical properties between the different dust 
species are likely to become apparent. This may explain 
at least partly the differences in distributions between the 
obscured M-type AGB stars and obscured carbon stars in 
the "Solar neighbourhood" . 

Obscured M-type AGB stars and obscured carbon 
stars seem to show very similar dependencies of their dust- 
to-gas ratios and especially their mass-loss rates on initial 
metallicity. There is no a-priori reason why the mass-loss 
mechanism should have a different dependency on initial 
metallicity for stars with different chemical types of cir- 
cumstellar dust, other than due to a different dependency 
of their dust-to-gas ratios on initial metallicity. The sim- 
ilarity between the dependencies of the dust-to-gas ra- 
tio on initial metallicity is surprising. For obscured M- 
type AGB stars one may expect that the fraction of met- 
als that are available for the formation of dust particles 
scales directly with the oxygen abundance in the photo- 
sphere, which scales at least approximately directly with 
the initial metallicity of the star. Hence a direct propor- 
tionality between dust-to-gas ratio and initial metallicity 
may not come as a surprise for obscured M-type AGB 
stars. For obscured carbon stars, however, the situation 
is very different: carbon stars only become carbon stars 
after 3 rd dredge-up has enhanced the carbon abundance 
in the photosphere from [C/O] < 1 to [C/O] > 1. For ob- 
scured carbon stars it is crucial to know the photospheric 
abundances of both carbon and oxygen, because the car- 
bon is locked into CO molecules until oxygen exhaustion 
and hence only the carbon excess is available for dust for- 
mation. It was thought that at lower initial metallicity 
the lower oxygen abundance would make it easier for 3 rd 
dredge-up to raise [C/O] above unity, but the fact that 
no optically bright luminous carbon stars were found in 
the MCs meant that it is not that simple (Iben 1981). Not 
only is it poorly understood how 3 rd dredge-up depends on 
initial metallicity (but see Marigo et al. 1999), there is also 
a second important phenomenon active: carbon star for- 
mation is avoided as long as the stellar mantle is massive 
enough to yield pressures and temperatures at the bottom 
of the convective layer sufficiently high for processing of 
carbon into oxygen and nitrogen to occur (Hot Bottom 
Burning; Iben & Renzini 1983; Wood et al. 1983). Thus 
it remains to be seen how the carbon excess for obscured 
carbon stars depends on initial metallicity (see also van 
Loon et al. 1999a). The data and analysis presented here 
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suggest that the carbon excess for obscured carbon stars 
may be directly proportional to the initial metallicity. 

6. Summary 

In conclusion, the comparison between the dust-to-gas ra- 
tios and mass-loss rates of obscured AGB stars in the 
SMC, LMC, and Milky Way suggests that the dust-to- 
gas ratios in the outflows of obscured AGB stars depend 
approximately linearly on the initial metallicity This was 
suspected by Habing et al. (1994) but still awaited ob- 
servational support. The total (gas+dust) mass-loss rates 
of obscured AGB stars show a much weaker dependence 
on initial metallicity. No appreciable differences are found 
between these relations for obscured M-type AGB stars 
and obscured carbon stars. 
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Appendix A: New near-IR counterparts of IRAS 
sources in the LMC 

Periods of weather conditions that were too poor for long-term 
photometric monitoring at the South African Astronomical 
Observatory (SAAO) at Sutherland, South Africa, in Decem- 
ber 1997 were used to search for near-IR counterparts of IRAS 
point sources in the direction of the LMC. This was done on 
the 1.9 m telescope with the Mk ill scanning photometer in the 
K-band. An aperture of 12" was used, chopping and nodding 
with a throw of 30". The search was limited to objects brighter 
than K = 13 mag. The areas around five IRAS point sources 
suspected to be obscured AGB stars (paper I) were searched. 
The candidate near-IR counterparts found are listed in Ta- 
ble Al, where the photometry is in the SAAO system (Carter 
1990), i.e. the J-band magnitude is transformed to the 0.75 
m telescope system. One object was re-observed under good 
photometric conditions (LI-LMC1284), together with the star 
HR2015 (8 Dor) for photometric calibration. Positions have 
been estimated by comparing the position of the diaphragm in 
the (red) acquisition video images with the second generation 
Digital Sky Survey, and are accurate to ~ 2". 

I retrieved 12, 25 and 60 /im data from the IRAS data base 
server in Groningen^j (Assendorp et al. 1995). Point sources 
were recovered by means of 2 x 2 square degree maps with 0.5' 

1 The IRAS data base server of the Space Research Organisa- 
tion of the Netherlands (SRON) and the Dutch Expertise Cen- 



pixels. The flux density was measured from a trace through 
the position of the source using the command SCANAID in the 
Groningen GIPSY data analysis software. LI-LMC203 shows a 
hint of duplicity: two similarly bright sources separated by 
one arcminute. LI-LMC987 looks slightly extended, and LI- 
LMC1284 is on top of brighter emission. LI-LMC1522 and 
especially LI-LMC1795 are isolated. Assuming identification 
of near-IR and IRAS sources, the spectral energy distribution 
were integrated graphically to yield bolometric magnitudes. 

LI-LMC203 is not identified with certainty. The best spatial 
coincidence is for the first listed in Table Al, that has blue near- 
IR colours incompatible with mass-losing AGB stars. There 
are two much brighter near-IR sources with moderately red 
(J — K) nearby, of which the third listed in Table Al is a 
cluster of ~ 10 stars within a diameter of ~ 12". The proposed 
near-IR counterparts of LI-LMC987 and 1795 have near-IR 
colours consistent with red giants without mass loss and are 
probably not associated with the IRAS sources. The near-IR 
counterpart of LI-LMC1284 is a heavily obscured AGB star, 
with IR colours compatible with either carbon or oxygen-rich 
dust. LI-LMC1522 is also identified as a dust-enshrouded star, 
with IR colours suggesting oxygen-rich dust. 



Appendix B: Expansion velocities 

Here the expansion velocities of AGB stars are briefly discussed 
in order to arrive at a justified calibration of the v cxp (L) scaling 
relation (Eq. (3)) for LMC stars. I consider v cxp derived from 
the separation of the peaks of OH maser line profiles, and from 
the width of CO (1-0) emission for Milky Way stars without 
OH measurements. The latter are divided by 1.12, following 
Groenewegen et al. (1998; see also Lewis 1991). 

The foxp is plotted against P in Fig. Bl. Periods for the 
LMC stars are from Wood et al. (1992) and Wood (1998), and 
their w cxp are from Wood et al. (1992) and van Loon et al. 
(1998b). Carbon stars (filled symbols) are distinguished from 
oxygen-rich, M-type stars (open symbols). As an AGB star 
evolves, P increases. Mass loss becomes important for PJJ400 
d and reaches a maximum for PJJ600 d (cf. Jura 1986). This 
evolution is also reflected in v cxp , which increases when PJ;400 
d and levels off at a value v cxp ~ 18 km s _1 for PiJ600 d (dotted 
line in Fig. Bl, see also Lewis 1991). The individual data sets 
(and Lewis 1991) suggest that AGB stars with PJS350 d evolve 
at ^constant u cxp typically 3 to 4 km s . This is also seen in 
(intrinsic) S-type stars with semi-regular variability (Fig. 18 in 
Jorissen & Knapp 1998). This may be a slow wind supported 
by another mechanism, such as radiation pressure on molecules 
(cf. Sahai & Liechti 1995; Steffen et al. 1997, 1998). It agrees 
with the observation that SiO masers, that depend on shocks 
by stellar pulsation (Alcolea et al. 1990), may be present for 
PiJ350 d and become ubiquitous when P<;700 d (Izumiura et 
al. 1994). 

Carbon stars appear to have somewhat larger v cxp than 
M-type stars at the same P, although the (few) SGC carbon 
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Table 3. Near-IR stars near IRAS point sources in the direction of the LMC (LI=LI-LMC: Schwering & Israel 1990) 
that are candidate obscured AGB stars. Listed are IRAS flux densities (in Jy), near-IR position, distance to the IRAS 
source (in arcsec), near-IR magnitudes, and bolometric magnitude assuming association (see text). Values between 
parentheses are l-a errors. 



LI 


Fl2 


F25 


F 60 a(2000) 


(5(2000) 


A 


J±ctj 


H±a H 


K±a K 


L±a L 


M hol 


203 


0.31 


0.20 


4 55 40.6 


-69 26 40 


8 


11.955(0.017) 




11.858(0.039) 




-5.8 








4 55 41.7 


-69 26 22 


26 


8.776(0.016) 




7.719(0.021) 




-8.3 








4 55 35.7 


-69 26 56 


22 


8.455(0.017) 




7.346(0.022) 


6.850(0.073) 


-8.6 


987 


0.48 


0.43 


5 24 41.7 


-69 15 20 


15 


13.006(0.069) 


11.906(0.035) 


11.622(0.053) 




-5.7 


1284 


0.37 


0.79 


5 32 38.9 


-68 25 22 


13 




14.015(0.049) 


11.896(0.023) 


9.127(0.071) 


-5.8 


1522 


1.00 


0.92 


5 40 13.1 


-69 22 50 


3 


10.660(0.015) 


9.919(0.018) 


8.612(0.017) 


6.886(0.032) 


-7.4 


1795 


0.25 


0.28 


0.16 5 56 42.7 


-67 53 21 


29 


12.058(0.025) 


11.157(0.019) 


10.907(0.024) 


10.336(0.129) 


-5.7 



stars do not obey this trend. The data is also consistent with 
smaller P for carbon stars at the same w cxp . 

SGC stars within 1 kpc from the galactic plane (squares) 
are distinguished from SGC stars beyond that (large circles). 
The latter are presumably of sub-solar initial metallicity and 
have smaller u cxp than the former. 

Wood et al. (1992) showed that « cxp is smaller at LMC 
metallicity than at solar metallicity, providing supportive ev- 
idence for Eq. (3). The LMC star with w oxp = 24 km s _1 
is IRAS04553— 6825, a very luminous RSG (van Loon et al. 
1998b, and references therein). The other LMC stars are AGB 
stars with « cxp ~ 10 km s _1 . The OH/IR stars in the Groe- 
newegen sample with PJ;1000 have u oxp ~ 10 to 20 km s~ , 
suggesting that initial metallicities of these Milky Way stars 
are higher than the LMC stars. The expansion velocities of 
the obscured AGB stars in the Galactic Centre range up to 
w cxp = 25 km s _1 , and initial metallicities of two to three times 
solar have been suggested by Wood et al. (1998). Considering 
all this, Eq. (3) is calibrated by demanding a star with LMC 
metallicity and L — 30, 000 L Q (M bo i = —6.5 mag) to have 
u cxp = 10 km s _1 . 
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Fig. B.l. Expansion velocity versus pulsation period, for samples in the LMC (Wood et al. 1992; Wood 1998; van 
Loon et al. 1998b), South Galactic Cap (Whitelock et al. 1994), Galactic Centre (Wood et al. 1998) and "Solar 
neighbourhood" (Groenewegen et al. 1998). Carbon stars are represented by filled symbols. The dotted line is drawn 
to guide the eye: expansion velocities are low for P<350, increase for longer periods, and reach a more constant level 
for P<;600. Low-metallicity stars in the LMC and SGC samples (at more than 1 kpc from the galactic plane) have 
lower w exp than stars with ^solar metallicity. 
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